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in : These lectures discuss how the direct detection of gravitational waves can be used to 

probe the very early Universe. We review the main cosmological mechanisms which 
could have produced relic gravitational waves, and compare theoretical predictions 
with capabilities and time scales of current and upcoming experiments. 

m : 
o ■ 

— ■ 1. Overview of gravitational-wave research 

In 1916 Einstein realized the propagation effects at finite velocity in the 
5— i , gravitational equations and predicted the existence of wave-like solutions 

of the linearized vacuum field equations P^. The work of Bondi |2] in the 
mid 50s, applied to self-gravitating systems like binaries made of neutron 
stars and/or black holes, proved that gravitational waves (GWs) carry off 
' energy from those systems, as the 1974 discovery of the binary pulsar PSR 

1913+16 by Hulse and Taylor |3] confirmed indisputably. This discovery 
is an indirect observation of GWs. In fact, the speed up of the two-body 
orbital period can be explained as due to GW emission. 

The experimental search for direct observation of gravitational waves 
begun only in the 60s with the pioneering work of Joseph Weber, and for 
almost three decades it has been pursued solely using meter-scale resonant- 
bar detectors. These detectors are cylindrical-shape bars whose mechanical 
oscillations can be driven by GWs. Five cryogenic resonant-bar detectors 
are in operation since 1990: ALLEGRO, AURIGA, EXPLORER, NAU- 
TILUS and NIOBE 0]. They are narrow-band detectors sensitive to GW 
frequency ~kHz. During the last years a world-wide network of kilometer- 
scale ground-based laser-interferometer detectors has been built and has 
begun operation in Japan (TAMA 300), in the United States (LIGOs) and 
Europe (GEO 600 and VIRGO) The frequency band is ~ 1-10 3 Hz. 
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Meantime the first detectors begin the search, development of the next gen- 
eration kilometer-scale interferometers is already underway |H1 El E3j • 
A laser-interferometer space antenna (LISA) [111 Ej is also planned by the 
European Space Agency (ESA) and NASA, though not yet fully funded, 
and could fly in ~ 2011. It consists of three drag-free spacecrafts in a 
triangle configuration. The spacecraft tracks the distance to each other (~ 
5 million kms) using laser beams, searching GWs in the frequency range 
~ 10~ 4 1CP 1 Hz. Thanks to the theoretical and experimental progress 
made by the GW community during the last forty years, the direct detec- 
tion of GWs is not far ahead of us and, hopefully, it will mark the next 
decade. 

The most promising sources for both ground- and space-based detectors 
are astrophysical sources at relatively low red-shift > like binaries made 
of black holes (BHs) and/or neutron stars (NSs), white-dwarf (WD) bi- 
naries, supermassive black holes, low-mass X-ray binaries, supernovae and 
pulsars. The large scale interferometers were indeed conceived and designed 
to detect GWs from those sources. These lectures will focus on cosmolog- 
ical sources at much higher red-shift z > 1. We shall discuss if and how 
the detection of GWs can be used to investigate physical processes occured 
around Universe's birth. Various excellent lectures and reviews were writ- 
ten on this subject like, e.g., the ones by Allen J3| an d Maggiore ^3] [ see 
also Sections 2.9 and 3.6 by Cutler and Thorne |12|], 

It is not possible to cover in two lectures all the literature on the subject, 
so we have to make a selection. In the first lecture, after briefly reviewing 
the key ideas underlying GW detectors [see Section 0, we discuss general 
features of relic GW spectra, their typical intensity and range of frequencies 
and phenomenological bounds [see Section EJ. The bulk of the lecture will 
be the analysis of the phenomenon of amplification of quantum-vacuum 
fluctuations, and the evaluation of stochastic GW backgrounds in stan- 
dard, quintessential and superstring-motivated inflationary scenarios [see 
Section 0] . The second lecture will review, in Section other physical 
mechanisms which could produce GWs, like first-order phase transitions, 
turbulence, preheating, etc., and cosmic strings [see Section Ej. Then, in 
Section we discuss which differences in the relic GW spectrum we could 
expect from brane-world scenarios and, finally, in Section|S|we comment on 
the possibility of extracting the red-shift - luminosity curve and cosmological 
parameters by detecting GWs from binary systems. Section El summarizes 
the conclusions. 
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2. Key ideas underlying GW detectors 

The simplest GW detector we can imagine is a body of mass m at a distance 
L from a fiducial laboratory point, connected to it by a spring of resonant 
frequency coq and quality factor Q. Einstein equation of geodesic deviation 
predicts that the infinitesimal displacement AL of the mass along the line 
of separation from its equilibrium position satisfies the equation |15j (valid 
for GW wavelengths 3> L and in local Lorentz frame of the observer at the 
fiducial laboratory point) 



AL(t) + 2-^ AL(t) 



ujI AL(t) = £ F+ h+(t) + F x h x (t) 



(1) 



where -F+, x are coefficients of order unity which depend on the direction of 
the source and the GW polarization angle; h + ^ x are the two independent 
polarizations of the GW. 

Laser-interferometer GW detectors, such as GEO, LIGO, VIRGO and 
TAMA, are composed of two perpendicular kilometer-scale arm cavities 
with two test-mass mirrors hung by wires at the end of each cavity. The tiny 
displacements AL of the mirrors induced by a passing- by GW are monitored 
with very high accuracy by measuring the relative optical phase between 
the light paths in each interferometer arm. The mirrors are pendula with 
quality factor Q quite high and resonant frequency luq much lower (~ 1 Hz) 
than the typical GW frequency (~ 100 Hz). In this case Eq. QJ, written in 
Fourier domain, reduces to AL/L ~ h. The typical amplitude, at 100 Hz, of 
GWs emitted by binary systems in the VIRGO cluster of galaxies (15 Mpc 
distant), which is the largest distance the first-generation of ground-based 
interferometers can probe, is ~ 10 -21 . This means AL ~ 10~ 18 m, a very 
tiny number! It seems rather discouraging, especially if we think to monitor 
the test-mass motion with light of wavelength nearly 10 12 times larger. It 
took a long time, theoretically and experimentally to develop sophisticated 
technology to achieve measurements of such tiny displacements. 

The electromagnetic signal coming out the interferometer will contain 
the GW signal but also noise - for example thermal noise in the suspension 
system and in the mirror itself, can shake the mirror mimicking the effect of 
a GW. The root-mean-square of the noise is generally expressed in terms of 
the noise power per unit frequency S n by the relation h lma ~ y S n (f) Af ~ 
AL/L, AL being the mirror displacement induced by noise and Af the 
frequency bandwidth. 
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3. Production of relic gravitational waves: general features 
3.1. Gravitational-wave spectrum 

Since in the following we will focus mainly on the primordial stochastic 
background of GWs, in this section we relate its energy density to the 
spectral density Sh, which is the quantity of direct interest when we want to 
compare theoretical predictions to experimental sensitivities. [Henceforth, 
we pose H = 1 = c] 

The intensity of a stochastic background of GWs at (present) frequency 
/ is generally expressed in terms of 

M/)-^, (2) 

Pc dlogf 

where p c = 3H 2 / (8itGn) and paw is the energy density of the GW stochas- 
tic background. a In the transverse-traceless (TT) gauge, denoting with ft 
the direction of propagation of the GW, with to, n the unit vectors form- 
ing with ft an orthonormal basis, and with = rhirhj — nihj 7 efj = 
TOj fij +hi rhj the polarization tensors, i, j = 1,2,3 [e^(^) e Qlj (£!) = 2(5 PQ ], 
a stochastic, isotropic GW can be written as 

/+oo p 
df / dtlhp&tye-W* (3) 
r = + ,x -°° J 

where dft = d<fidcos8, ft,p(/, ft) = hp(— /, ft), and in writing Eq. © wo 
have assumed, in full generality, that the wave is located at x = 0. If 
the stochastic background is also stationary and unpolarized, averaging the 
Fourier components over an ensemble gives |14| 

< h P (f, ft) h Q (f, ft') >= 6(f - /') s 2 (n, ft') i- s PQ 1 s h (f) , (4) 

with <5 2 (ft, ft') = 8(cos9 — cos 6') 8{cj) — (/)'). The above equation defines the 
(one-sided) spectral density Sh [Sh(f) = Sh{~f)]i which has the dimensions 
of Hz -1 . The energy density of GWs is given by the 00 component of 
the energy-momentum tensor averaged over several wavelengths [see e.g., 
Ref. [SI] 



a The value of Ho is generally written as Ho = ho X 100 km/sec/Mpc where ho param- 
eterizes the existing experimental uncertainty. Henceforth, all primordial spectrum will 
be expressed in terms of h^ SIgw- 
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Figure 1. Sensitivities, expressed in terms of SIqWi versus frequency of LISA and 
ground-based detectors of first, second and third generation. On the top axis, the es- 
timated temperature of the Universe when GWs are produced by causal mechanisms 
during RD era [see Eq. I1UI 1 , 



Inserting in the above equation Eq. @ and using the property that for a 
stochastic background the spatial average can be traded with the ensemble 
average defined by Eq. (@J, a straightforward calculation leads to: 



4?r 2 1 / f \ 2 

"gw(/) = xt72/ 3 ^(/)^ 1.25 x 10 36 -2M- fS h (f). (6) 



o 



Note that if ilcw ~ C f a with < a < 3, due to the factor f 3 appearing 
in the RHS of Eq. © , a GW detector operating at lower frequency needs 
lower sensitivity than a GW detector operating at higher frequency. This 
explains why space-based detectors like LISA can have better performances 
than ground-based detectors, although the expected sensitivity of LISA is 
worse than LIGO and VIRGO. 

In Fig. H w e plot the sensitivities for ground-based detectors of first 0, 
second (~ 2008) [3 IHJ and third (> 2011) HHHj generation. More specif- 
ically, for the second generation we use the Advanced LIGO configuration 
or LIGOII [ZIIHj; for which considerable research and development (R&D) 
have been gained during the last years, and for the third generation we 
use, as an example, the speed- meter configuration The R&D for 

third generation of interferometers as LIGOIII/EURO started only more 
recently. 
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A stochastic background is a random process which is intrinsically in- 
distinguishable from the detector noise. As we shall see in the following, 
the GW signal is expected to be far too low to exceed the noise level in any 
existing or planned single detector on the earth. Moreover, the instrumen- 
tal noise level will not be known sufficiently well a priori to search for excess 
noise in each instrument. Therefore, the optimal strategy which has been 
proposed is to perform a correlation between two or more detectors, 
possibly widely separated to minimize common noise sources. By correlat- 
ing two detectors the increase in sensitivity h lrns is ~ (A/T) 1 / 4 where A/ 
is the bandwidth and T the total observation time. Henceforth, we shall 
show the sensitivities obtained correlating two ground-based detectors. As- 
suming a constant GW spectrum, and correlating for fourth months, we 
obtain at 90% confident level for two LIGOI: /ig £!gw — 3.5 x 10~ 6 , for two 
LIGOII: h% n GW ~ 5.1 X 1(T 9 and for two LIGOIII: h% fl GW ~ 3.7 x 10~ u . 

In the case of LISA, since only one space-based detector is currently 
planned, the method of correlating two detectors cannot be used. [In 
Ref. |17| optimal orbital alignements for correlating a pair of future 
LISA-kind detectors have been investigated.] Armstrong, Estabrook and 
Tinto JSl realized that by combining the signals from the three spacecrafts 
which compose LISA, it is possible to measure the instrumental noise power. 
Then, any excess noise above this power will be due to GWs. However, 
besides the primordial stochastic GW background, there are various as- 
trophysically generated stochastic GW background due to binaries present 
in our galaxy or outside it. These stochastic backgrounds are produced 
when the incoherent superposition of gravitational radiation emitted by a 
large number of astrophysical sources cannot be resolved individually. A 
possible, although rather challenging way of discriminating between the 
primordial and astrophysically generated backgrounds was suggested by 
Ungarelli and Vecchio ^Hj- As LISA rotates along its orbit, its sensitivity 
to different directions changes, so LISA could measure the isotropy of the 
stochastic background and use it to separate the primordial stochastic back- 
ground, which is isotropic, from the galactic stochastic background, which 
is anisotropic, being concentrated mostly in the galactic plane (rather than 
in the halo). In Fig. ^ we plot the planned LISA sensitivity 40 . 

3.2. When gravitons decoupled: thermal spectrum? 

It is well known that particles which decoupled from primordial plasma at 
time tdec, when the Universe had temperature Td oc , carry information on 
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the state of the Universe at T^ cc . The weaker the interaction, the higher the 
energy scale when they drop out of thermal equilibrium. To estimate Ta ec 
for gravitons, we proceed as follows [201 ED ■ We assume that gravitons are 
in thermal equilibrium in the primordial plasma through point-like four- 
body interactions involving two gravitons. We denote with Y = na\v\ 
the interaction rate per particle, n being the number density of particles 
in equilibrium, a the scattering cross section, and v the relative velocity. 
From dimensional considerations we expect a ~ G^ E 2 ~ G^T 2 , where 
Gn is the Newton constant and E 2 is the average energy squared. For 
relativistic particles in equilibrium at temperature T, n ~ T 3 . We also 
assume v ~ 1. The interaction rate is then Y ~ G^T 5 . Assuming the 
Universe evolves adiabatically, we have gsT 3 a 3 = const., where gs is the 
number of effectively massless degrees of freedom [sec Eq. (3.73) of Ref. 21 ]. 
Since for relativistic particles p ~ T 4 , we have T/T oc H oc T 2 /Mpi, thus 
the condition to maintain thermal equilibrium is Y > i7, i.e. the interaction 
time-scale must exceed the local Hubble time. This gives: 




So, as it could have been anticipated from dimensional considerations, gravi- 
tons decoupled at Td cc ~ Mp\. Let us derive at which temperature it cor- 
responds today. The temperature of the Universe, i.e. of the particles in 
thermal equilibrium, scales as T oc g s a -1 , while the temperature as- 
sociated to particles which have decoupled from the thermal bath drops 
as T oc a -1 . Applying these considerations to the thermal bath of gravi- 
tons, and using the fact that at T > 300 GeV the Standard Model (SM) 
of particle physics predicts gs — 106.75, while at present time, assuming 
three neutrino species, gs — 3.91, it is easily derived that the graviton tem- 
perature today is at most [Hj T g = (3.91/106. 75) 1 / 3 2.75 K ~ 0.9 K [at 
most because the effective number of massless degrees of freedom could be 
higher then 106.75 at Planckian energies.]. Thus, the thermal spectrum of 
gravitons ranges in the GHz and is given by |13| 

"•<" (i)' i 37#TT — >h °' (i)' ■ <8 > 

with f g = k T g /h ~ 1.9 x 10 10 Hz. [Note that for CMBR we have T 7 ~ 2.73 
K and / 7 = 5.7 x 10 10 Hz.] However, since gravitons interact very weakly 
with matter and radiation, the time required to reach thermal equilibrium 
could be longer than the Hubble time. So, it is quite unlikely that this 
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radiation ever existed. Moreover, if the Universe underwent a period of 
inflation after the Planck era, the background should have been strongly 
diluted. More importantly, the theory of gravity at the Planck energy could 
differ significantly from Einstein theory; if so, the above analysis should be 
modified, accordingly. 

3.3. Gravitational waves produced by casual mechanisms: 
typical frequencies 

Two features determine the typical frequency of GWs of cosmological origin 
produced by some casual mechanism |14j : (i) the dynamics, which is model 
dependent, and (ii) the kinematics, that is the red-shift from the production 
era. 

Let us assume that a graviton is produced with frequency /* at time 
t* during RD or MD era. We want to evaluate its frequency today 
/ = /* a*/ao- Assuming that the Universe evolved adiabatically, that is 
9s(T*)T? al = g s (T ) T§ a 3 , and using T = 2.73K and g s = 3.91 we get 



Since the size of the Hubble radius is the length scale beyond which causal 
microphysics cannot operate, from causality considerations we expect the 
characteristic wavelength of gravitons produced at time t* is either ~ H~ x 



or smaller. So, we pose |2] A* = eH^ 1 with e<l. If the GW signal is 
produced during RD era H% = 8ttG n Prad/3 = 87r 3 ff * Tf /(90 Mpj) and 



On the top axis of Fig.^we show the production temperature T» as obtained 
from Eq. (|lt)(l by expressing T* in terms of the frequency /. For the kind of 
mechanism analysed in this section, LISA could probe physics at TeV scale, 
ground-based interferometers in the range ~ 10 8 -10 10 GeV, while GUT 
and Planck scales could be probed in GHz region, where however no GW 
detectors have been planned so far, but electromagnetic (EM) microwave 
cavities have been proposed 

3.4. Phenomenological bounds 
3.4.1. Big-bang nucleosynthesis bound 

The theory of Big-bang nucleosynthesis (BBN) j2H] predicts rather success- 
fully the primordial abundances of light elements like deuterium, 3 He, He 
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(10) 
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Figure 2. Summary of phenomenological bounds on energy density of relic GWs. 



and 7 Li. If at nucleosynthesis time (T ~ MeV), the contribution of the pri- 
mordial GWs to the total energy density is too large, then the expansion 
rate of the Universe H , and the freeze-out temperature which determines 
the relative abundance of neutrons and protons, will be too high. Thus, 
neutrons will be more available and 4 He will be overproduced, spoiling 
BBN predictions. Detailed calculations provides the following bound on 
the energy density in GWs integrated over frequency, 

/•/=+oo 

/ dlogfh 2 n GW (f) <5.6x l(T 6 (A^-3), (11) 

Jf=o 

where N u is the effective number of neutrino species. The bound (flip holds 
only for GWs that were already produced at time of nucleosynthesis. Since 
the integral cannot exceed the bound, also its positive definite integrand 
cannot do it, unless flow <C 10~ 6 for most of the frequency range and 
has a very narrow peak on the order of ~ 10~ 5 at some frequency. This 
possibility is however not very plausible. Thus, Eq. (|llfl gives for A/ ~ / 
and N„ < 4 the value hi fl GW < 5.6 x 10~ 6 , which is plotted in Fig. El 



3.4.2. COBE bound 

The COBE bound comes from the measurement of temperature fluctua- 
tions in the Cosmic Microwave Background Radiation (CMBR). Indeed, 
a background of GWs at very low frequencies can produce, if sufficiently 
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strong, a stochastic red-shift of the CMBR frequencies, through the well 
known Sachs- Wolfe effect. Indicating by ST the fluctuation induced in the 
CMBR temperature, we have 



"gwCO < -r Hf 3 x 10-18 Hz < f < 10-16 Hz > ( 12 ) 



where the range of frequencies is determined by the condition the GWs 
are inside the Hubble radius today (/ > Hq ~ 3 x 10~ 18 Hz) and 
they were outside the Hubble radius at the last scattering surface (LSS) 
(/ < y^hssH ~ 10 -16 Hz with z LSS ~ 10 3 ). The value of 5T/T induced 
by GWs on CMBR cannot exceed the observed value 5T/T ~ 5 x 1(T 6 . 
Detailed analyses give [2U Ej : 



hi n GW (/) < 7 x HT 11 3 x l(T 18 Hz < / < lCr 16 Hz. (13) 



GWs can saturate this bound if the contribution from scalar perturbations 
is subdominant, and this depends on the specific inflationary model. The 
COBE bound (S) is shown in Fig. 

3.4.3. msec pulsar bound 

The very accurate timing of msec pulsars constraints Hqw- in fact, if a 
GW passes between us and the pulsar the time of arrival of the pulse is 
(doppler) shifted. Many years of observations lead to the bound |2H IT1) 



where /* is derived from the total observation time of ~ 8 years [see Fig.E]. 

4. Amplification of quantum vacuum fluctuations 

The amplification of quantum vacuum fluctuations is a very general mech- 
anism characterizing quantum field theory in curved space time, first dis- 
cussed in cosmology by Grishchuk and Starobinsky [2E]- Let us first 
give a formal description of it. 

We assume that the background field dynamics is described by the ac- 
tion: 






/ > /„ ee 4.4 x KT 9 Hz , 



(14) 




(15) 
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where 1Z is the Ricci scalar, S m refers to possible matter sources with 
\/\g\T flv — 2SS m /5g' J ' v , T M „ being the energy-momentum tensor. We 
restrict ourselves to an isotropic and spatially homogeneous Friedmann- 
Lemaitre-Robertson- Walker (FLRW) unperturbed metric, with scale fac- 
tor a and ds 2 — g^ v dx^ dx v = —dt 2 + a 2 (t) dx 2 . The free linearized wave 
equation for the TT metric perturbations Sg^ = (h^o = 0, V,,/!^ = 
0, h^n = 0) can be obtained by perturbing Einstein equations and keeping 
all the sources fixed, 8T^ = 0. A straightforward calculation gives: 

ah?(t,x) = -L=d ll (y/\i\ ! rd v ) h/(t,x) =0. (16) 

Vm 



Introducing the conformal time 77, with dr\ = dt/a{t) and writing 

P=+,x £ 

each polarization mode h^(r]) satisfies the equation 

h>l( V ) + 2^hi~( v ) + k 2 h j :( v ) = 0, (18) 

where we indicate with a prime the derivative with respect to conformal 
time 77. 

For simplicity let us consider two cosmological phases I and II, e.g., in- 
flation and RD or MD phase. We denote with h^(r]) and H^(rj) the solution 
of Eq. l|18f) when the scale factor of phase I and II is used, respectively. The 
mode expansion of h? in phase I reads: 

hi = ^G^Y. J V2fc(2^ [ fo pWM^f j (f))e^ + H.C.] , (19) 

where H.C. stands for hermitian conjugate, b\>(k) is the annihilation opera- 
tor with respect to the vacuum in phase I, i.e. 6p(fc)|0)i = with P = x , + 
[classically, it can be traded with a random variable.] We refer to x and k 
as the comoving coordinate and momentum, related to the physical quan- 
tities by Xphys = a(t) x and fc p i iys = k/a(t). The mode expansion of in 
the phase II can be obtained from Eq. (fH^I with the substitution hj: — ► 
bp(k) — » bp(k). The annihilation operator of phase II satisfies the equations 
bp(k)\0) u = with P = x,+. The quantities h%} and {H^HZ} are 
complete bases and can be related to each other by the so-called Bogoliubov 
transformation 
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%' 

Note that the above equations express the positive-frequency components 
of the operator in phase I in terms of both the positive and negative com- 
ponents of the operators in phase II, and viceversa. The coefficients arp 
and (3?p are called Bogoliubov coefficients and satisfy the equations 



£ 

fc 

£ 

fc 



[^k^k ~ = S k,k 2 > ( 22 ) 

K^^-^fe^fefe] =°- (23) 



Inserting given by Eq. (|2(J[1 . in the expression of and equating 
to ft./, we find the well-known relation between annihilation and creation 
operators in the two phases: 

bl W = E [ a Pk bU (k')+Phb U Hk')] , (24) 
%' 

b l \k) = £ [a^ b\k') - % 6 I t(fc')] • (25) 

Thus, if /3gg, 7^ 0, even starting from the vacuum i(0|6 I " i '(fc) 6 I (fc)|0)i = 0, 
we end up with a final number of particles given by 

I (0|& II t(fc)& II (fc)|0) I =^| i 3^,|V0. (26) 
%' 

This phenomenon is known as amplification of quantum-vacuum fluctua- 
tions in curved space time. It is due to the inevitable presence of positive 
and negative components in the basis I when expressed in terms of the ba- 
sis II [see Eqs. (|20fl ]. The existence of the two complete bases, which are 
associated to two independent Fock spaces, is a consequence of being the 
space-time curved. 

The physical idea underlying the quantum-vacuum production of par- 
ticles can be also explained as follows [Henceforth, we assume an 
isotropic and spatially homogeneous metric, so we can pose a^g, = at 8zz, 
and fizz, = Pk5^%i] Suppose that over a time-scale AT ~ H^ 1 the cos- 
mological phase changes, e.g., from phase I to phase II. If f* is the time 
at which the transition occurs and /* is the physical frequency of the 
mode at that time, there exist two possible regimes: abrupt transition 
for which 2-irf* if" 1 <C 1 and adiabatic transition 27r/* ff" 1 ^> 1. Let us 
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assume that the quantum state \q) before the transition has number of 
particles Nj. = (<z|&i b\\q). For modes for which the transition is sudden, 
i.e. 27T/* H^ 1 <C 1 (for these modes A ^ H^ 1 , that is their wavelengths 
are larger than the Hubble radius or in jargon they are said to be outside 
the horizon, the physical state does not have time to change during the 
transition. However, after the transition, the number of particles should be 
evaluated from annihilation and creation operators of phase II, 

N? = (q\b^ bf\q) =Nl(l + 2|/? fc | 2 ) + |/3 fc | 2 . (27) 

Therefore, if the number of particles when inflation starts is Nl , this number 
is amplified by the factor 1 + 2|/3fc| 2 jSHllEl- Moreover, due to the mixing 
between positive and negative frequencies even the vacuum state of phase I 
(Nl = 0) is a multiparticle state when referred to the vacuum state of phase 
II. On the other hand, for modes for which the transition is adiabatic, i.e. 
271-/* H* 1 ^> 1 (for these modes A<C H^ 1 , so their wavelengths are smaller 
than the Hubble radius or in jargon these modes are said to be inside the 
horizon), the quantum state has the time to follow the evolution of the 
scale factor. No particle production occurs in this case. In some sense the 
modes for which the transition is adiabatic, do not feel any effect of being 
the spacetime curved or changing in time. No mixing between positive 
and negative frequencies occurs, as it is in flat spacetime. [Practically, 
the relic GWs spectrum drops off exponentially for frequencies larger than 
ycutoff ^ if*/(27r), i.e. for modes which never went outside the Hubble 
radius.] 

We want now to give a more intuitive explanation of the phenomenon 
of amplification of quantum- vacuum fluctuations, which is very close to the 
original derivation due to Grishchuk 23 an d Starobinsky It is based 
on a semiclassical approach. By introducing the variable ipkiv) — Q-hkiv)? 
hk being one of the two polarization modes, Eq. (|TB)l can be recast in the 
form 

i/> + [k 2 - U(n)} V fc =0, U(n) = -, (28) 

a 

which is the equation of an harmonic oscillator in the time-dependent po- 
tential U(rj). [Note that to get unambiguous results, the change of variables 
ijj(k) — > h(k) can be used only if the scale factor and its first derivative are 
continuous at the transition between phase I and II .] For simplicity, 
let us consider a de Sitter inflationary era, i.e. a = —l/(r)Hds), H^s being 
the constant Hubble parameter during de Sitter phase. 
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If k 2 < \U(r))\, since \U(if)\ ~ ?r 2 and (a#ds) ~ 1/ry, then fcr; < 1, 
which means fc/a <C -ffds or H^g f, i.e. the mode is outside the Hubble 
radius (or in jargon under the potential barrier |?7(?7)|). In this case the 
solution of Eq. (|2"%)l is 



A k + B k 



dt] 



a 2 (r)) 



A »- B > iJt)- (29) 



Since during de Sitter era the scale factor gets larger and larger, the second 
term in the RHS of Eq. (|29(l for h k becomes more and more negligible, 
thus the tensorial perturbation h k remains (almost) constant while outside 
the Hubble radius. On the other hand, if k 2 » we have kr\ 3> 1 

or A <C H^g , which means the the mode is inside the Hubble radius (or 
in jargon over the potential barrier | £/(??) |). The solution of Eq. (|28|l is 
a plane wave ip k ~ e ±lkv . Thus, in this case the tensorial perturbation 
hk ~ e ±lkr> /a decreases in time while inside the Hubble radius. Therefore, 
in de Sitter case, the longer the tensorial-perturbation mode remains outside 
the Hubble radius, the more it gets amplified. When the fluctuations re- 
enter the Hubble radius, the solution of Eq. Ij28(l is 

V>fe~afce- 4fc " + /? fc e 4fc ", (30) 



and contains both positive and negative modes. The coefficient (3k is the 
Bogoliubov coefficient and it can be determined imposing the continuity of 
ipk and ip' k all along the cosmological phases. 

Before ending this section we want to express the intensity of the 
stochastic GW background J2J in terms of the number of gravitons per 
cell of the phase space, which we indicate by rif with / = |fc|/(27r). In 
the case of GWs produced out of vacuum, n/ = A'j 1 , thus it is the cru- 
cial quantity to evaluate. For an isotropic stochastic GW background 
Pgw = 2 jd 3 k/(2n) 3 (kn k ), thus 

fW/H-ie^n// 4 . (31) 

Pc 

From the above expression it is straightforward to deduce that if the 
planned space- and ground- based GW experiments will detect a stochas- 
tic GW background whose spectrum satisfies the BBN phenomenological 
bound discussed in Section l3.4.1l i.e. J]qw < 10~ 6 , the occupation number 
rif » 1. This means the observed stochastic GW background is classical. 
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4.1. Standard inflation 

In this section we shall first derive the GW spectrum in de Sitter inflation, 
which although is an ideal case (because it does not lead naturally to the 
transition from inflation to RD era), it has the advantage of providing 
analytical solutions for the tensorial perturbations. Secondly, we discuss 
how the GW spectrum is modified in more realistic inflationary scenarios, 
like slow-roll inflation. 

4.1.1. de Sitter inflation 

The scale factor during de Sitter era is a(rj) = — l/(H(\sr]) with — oo < 
rj < ry*, while during the RD era we have a(n) = (n — 2n*)/(Has'i]*) with 
?7* < rj < 77 q, where r] eq is the time at which there is equality between 
radiation and matter era. [Henceforth, we refer to quantities evaluated at 
present time with the subscript 0.] 

Let us observe that during de Sitter phase any (classical) tensorial 
perturbation is de-amplified and their wavelength gets always stretched. 
Indeed, for a given comoving wavenumber k, it can be easily found [23 
that the de-amplification coefficient is A(f) = h(f,r]o)/h(fi n ,r]i n ) = 
(H d sa*/(H a )), where / = k/(2ira ), f in = fc/(27ra in ) and Af in = 
log(a*/a; n ) is the number of e- foldings. [To solve the horizon problem |32) 
the minimal amount of e-foldings is Af m in — log(-ffds o-*/(Hq &o)).] Thus, 
the longer the inflationary era, the smaller the coefficient A. Moreover, the 
stretching of the wavelength is A/Ai n = e.^ (ao/a*). 

In de Sitter case Eq. (fT%f can be solved exactly and the solutions are 
rather simple, b they read: 

h k {v) = \ (i-^) e ~ ikV > ~™<V<V* (32) 

h k (77) = - [a k e~ lk 1 + (3 k e lk «] , ^ < r, < i hq . (33) 
a 

[Here, we have assumed that the initial state of the Universe is the Bunch- 
Davis vacuum [33], so in Eq. l|2*7ll we could unambiguously set = 0. See 
Ref. [31] where this assumption is relaxed.] The Bogoliubov coefficients a k 
and f3k in Eq. I|33f) are obtained imposing the continuity of h k and h' k across 
the transition. A straightforward calculation jSH] gives 3k = l/(2fc 2 ry^), 
so the occupation number is N k — |/3fc| 2 - Using the following relations: 



b These solutions can be obtained easily rewriting Eq. 1181 in terms of the function u^, 
defined by = (au^)' /a 2 . 
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fc|?7*| = 2nfa Q \ri*\ = 2n f a /(a* H ds ) where a* = l/(iJ ds |?7* |), a /a* = 
(io/t oq ) 2/3 (i cq A*) 1/2 , we obtain k\r^\ = f/f, with 

where we used z cq ~ 10 4 , i C q ~ 10 10 s and t* = The frequency /* 

is the cutoff frequency beyond which the GW spectrum falls down expo- 
nentially. By using Eq. 1|31[) with n& = l/(4fc 4 ?7 4 ), it is straightforward to 
find that the GW spectrum for fluctuations that re-enter the Hubble radius 
during the RD era is then 

^ GW (/),4 X 10-"( 6xl ^ M J 2 . (35) 

For fluctuations that re-enter the Hubble radius during the matter era, a 
similar calculation gives |33| : 

ftJnow(/)=!4x io- u {e^^Wif • < 36 > 

where / eq = H e J(2x)/(l + z cq ) ~ 10" 16 Hz. The COBE bound discussed 
in Section l3~4"2l imposes that the Hubble parameter during the de Sitter 
phase should be < 6 x 10" 5 M Pi |55| . 

4.1.2. Slow-roll inflation 

As a more realistic example we now discuss slow- roll inflation |32| . In this 
scenario a scalar field (the inflaton field) drives a period of accelerated 
expansion by rolling toward the minimum of its potential V ((/>). [Inflation 
can occur if the slow-roll conditions are satisfied 4> 2 <C V(4>) and \<fi\ <C 
\3H tj>\, \V'\ = \dV/d<f>\.] Differently from de Sitter inflation, in this case the 
Hubble parameter is not exactly constant during the inflation era, it slowly 
decreases in time, and as a consequence, the GW spectrum for fluctuations 
that go out of the Hubble radius during the inflation era and re-enter during 
RD era, is not flat but acquires a tilt: fto^Gw(/) ~ V /™ T > c where we 
denote by V the value of the inflaton potential at the time the present 
Hubble scale crossed the Hubble radius during inflation. The spectral slope 



c Whereas in de Sitter inflation the amplitude of the GW spectrum depends on the 
square of the (constant) Hubble parameter H^s [see Eqs. 1351 . 1361 1. in potential-driven 
inflation, it depends on V. Indeed, by applying the slow-roll conditions to the FLRW 
equations it can be easily seen that H 2 ~ 87tGn V/3. 
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Figure 3. In the left panel we sketch the evolution of the expansion rate of the Universe 
H and of some physical wavelengths during a period of standard inflation, followed by 
RD and MD eras. In the right panel we plot the stochastic GW background for slow-roll 
inflation obtained in Ref. 1361 for T/S = 0.18 and dn^/dlogk = (continuous line), 
T/S = 0.18 and dn T /d\ogk = -10" 3 (dot line), T/S = 0.0018 and dn T /d\ogk = 
(dash line) and T/S = 0.0018 and dnx/dlogk = — 10 -3 (dot-dash line) The sensitivity 
of space-, (correlated) ground-based detectors and the stochastic background from WD 
binaries is also shown for comparison. 



is riT = — (v' /V) 2 Mp[/(87r) and because of the slow-roll condition 
\tit\ <C 1. Quite nicely, the spectral slope can be expressed in terms of 
the scalar contribution S and tensorial contribution T = 0.61\^/Mpi to the 
quadrupole CMB anisotropy [see, e.g., Eqs. (3), (4) in Ref. 36 ]. The result 
is n T = -T/(7S). 

The explicit expression of the GW spectrum is given by Eqs. (5), (6) 
in Ref. In Fig. [21 we plot the spectrum for two values of T/S, taking 
also into account the first-correction for the variation of the spectral slope 
with scale, nr(fc), and compare it to the capabilities of LISA and (corre- 
lated d ) ground-based detectors. Figure shows that if the predictions of 
standard inflation are correct, unfortunately there is no hope of observing 
the stochastic GW background from slow-roll inflation with current and 
planned GW detectors. 

As seen from Fig. |3 in the case of LISA, primordial GW backgrounds 



d The curves have been produced by multypling LIGOl, LIGOII and LIGOIII sensitivities 
curves shown in Fig. Q by a constant factor obtained correlating Hanford (WA) and 
Louisiana (LO) LIGOs for four months at 90% confident level assuming a flat GW 
spectrum. 
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can be "covered" by the galactic WD binary background [37]. Even if two 
LISAs will be used, by correlating their data stream, the minimum achiev- 
able sensitivity is ft GW > lCT 13 37. Thus, to probe fl GW ~ 10~ 16 -10- 15 
from inflation with space-detectors of LISA-kind, the detectors have to op- 
erate outside the regime of mHz which is dominated by close WD binaries. 
Various scientists have already started discussing a follow-on mission of 
LISA. Given the current astrophysical understanding, the most promising 
region should be ~ 0.1-1 Hz 1381 lit)| HT] . which can be obtained by 
shortening LISA arms by a factor of 100. In this higher frequency band 
the astrophysically generated background from unresolved NS binaries is 
present; however, because the number of sources per frequency bin is ex- 
pected to be not very high, it should be possible to remove it. To achieve 
sensitivities on the order of JIqw ~ 10 _16 -10~ 15 , major, rather challenging 
improvements should be made, as increasing the laser power, the dimen- 
sions of the mirrors, etc.. 

The detection of relic GWs from inflation could be achieved also by 
measuring the polarization 02] m the CMBR with future CMB probes, as 
originally suggested in Ref. . Indeed, the polarization tensor on the ce- 
lestial sphere can be decomposed into its "gradient" and "curl" components. 
The scalar fluctuations do not have handness, so they do not contribute to 
the curl, while tensorial fluctuations do contribute to the curl. If the po- 
larization measurement will give a curl different from zero, GW will be 
detected and it will be possible to discriminate between the plethora of in- 
flationary models currently available. Indeed, some of those models predict 
that the tensorial component T is negligible. 

In Sections 14.21 14.31 we shall discuss other inflationary models. Before 
doing it, we want to point out that quite generally the tensorial perturba- 
tions will satisfy an equation similar to Eq. I|28|> but with the scale factor 
a replaced by the so-called "pump" field a. The pump field depends on 
the dynamics of the scale factor and the other fields which are part of the 
background. Following Refs. [441 1451 l5~4l 130) . if we have two cosmological 
phases with a = (?7/?y*) 71 for — oo < r/ < r/i < and a = \(ri — 2r/i)/r/i| 72 
for 7/i < 7] < r/2, then the solutions of Eq. 1)28(1 (having replaced a with a) 
can be expressed in terms of Hankel functions of first and second species, 
as 



Mri) = VW^^\[a k H£\k\r 1 -2 m \)+(3 k Hi 2 i \k\r 1 -2 m \)}, (38) 




(37) 
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where V\ — — 1/2| and v 2 = I72 — 1/2) , and we have normalized Eq. (|3*7I) 
allowing only positive frequencies at r\ — > —00, so that — ► Ce~ lkr> j \fk. It 
can be derived that the Bogoliubov coefhcient (3k entering the GW spectrum 
is given by [HOj 

l/3/| 2 ~/ 2£T er = 1- |7i| - I72I for 7i,72 > 1/2 or 7l , 72 <-l/2, 

(-T = — 72 1 all other cases . (39) 

The spectral slope parameter fij- is defined by 

dlogf2 GW . , . 

it = —r, 7— = 4 + 2e T , (40) 

"log/ 

and it is directly related to the kinematic behaviour of the background, 
that is to the evolution of the curvature scale and Hubble parameter. If 
we apply the above equations to de Sitter inflation, we find that 71 = — 1, 
for RD era 72 = 1 and for MD era 72 = 2. Thus, the spectral slope for 
fluctuations which re-enter the Hubble radius during RD era is tit = 0, and 
for those which re-enter the Hubble radius during MD era is = —2, in 
agreement with Eqs. (|35|) . 1136(1 . Moreover, in the case of slow-roll or power- 
law inflation, 71 < —1, and it can be easily found that for fluctuations that 
re-enter the Hubble radius during RD era tit = 2 + 2 7 i < 0, as anticipated 
in Section KHH 

Finally, we note that this way of estimating the spectral slope can be 
easily generalized to the case of n cosmological phases occurring at 
times rjx, t]2, • • • if the pump fields follow a power-law evolution with co- 
efficients 71,72,.... In this case the relic GW spectrum is characterized 
by various frequency regions or branches, and in each of them the spectral 
slope depends only on the kinematics of the phases in which the mode as- 
sociated to a given frequency went out of the Hubble radius and re-entered 
the Hubble radius. On the other hand, the amplitude of the GW spectrum 
and special features in it, like possible oscillations, are determined by the 
dynamics and by physical assumptions, such as when and how the reheat- 
ing process took place, whether entropy production is produced at some 
later stages, etc.. 



4.2. Superstring-motivated cosmology 

Potential-driven inflation of the kind discussed in Section 14.11 cannot be 
easily implemented in string theory if the dilaton field or a modulus field 
is simply identified with the inflaton field ■ This result forced people to 
conceive new ways of reconciling inflation and string theory. Henceforth, we 
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shall briefly discuss some of those attempts which use supergravity descrip- 
tion of superstring theory, the so-called pre-big-bang (PBB) scenario |4"T] 
and bouncing-Universe scenarios |48l I49j . Quite generally those scenarios 
predict a stochastic GW background anything but flat, which could be of 
interest for space- and ground-based detectors. 



4.2.1. Pre-big-bang scenario 

Due to the presence of other fields, at low energy superstring theory does 
not give Einstein general relativity - for example heterotic string theory in 
four dimensions is described by the action 



K + gT M V du<p - ^(dB) 2 - V(tp) 



(41) 



where ip is the dilaton field, related to the string coupling by g 2 = e v ; 
dB = d^Bjyp + d u B plJ/ + dpB^y, where i? M „ is the two-form gauge field or 
antisymmetric field; V(<fi) is a non-perturbative potential; and where A s is 
the string scale. In writing Eq. Q41JI we have disregarded for simplicity the 
internal dimensions, whose dynamics can be described in terms of moduli 
fields 02]. 

Hencerforth, we limit the discussion to the homogeneous and isotropic 
case with B = and V = [ds 2 = -alt 2 + a 2 (t) dx 2 ,(p = (p(t)}. In this case 
the solution of the low-energy string-effective action l|41|) satisfies the scale- 
factor-duality symmetry: a(t) — > l/a(t), tp(t) — > tp(t) — 61oga(t), c with 
a(t) ~ t 1 /^ anc [ n _ \ogt. Noticing this property, Veneziano |5U] 
conceived the idea of implementing the inflationary phase at times before 
the would-be big-bang singularity. Indeed, it is easily shown that for t < 0, 
a > 0,a > 0, thus the Universe undergoes a (super) inflationary phase. 
Two different but physically equivalent descriptions of the PBB phase exist: 
either the string-frame picture described by Eq. i|41|) . where the Universe 
undergoes an accelerated expansion (H > 0, H > 0, (f > 0), or the Einstein- 
frame picture, where the action (|41|l has the standard Hilbcrt-Einstein form 
and the evolution of the Universe is described by an accelerated contraction, 
or gravitational collapse (H < 0, H < 0, > 0). 

This new kind of inflation, which can be shown solves the homogeneity 
and flatness conundra, is driven by the kinetic energy of the dilaton field 
and forces both the string coupling (y > 0) and the spacetime curvature 
to grow toward the future. As a consequence, at least in the homogeneous 



e Here for convenience the origin of time has been fixed at t = 0. 
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case, the inflationary stage lasts for ever (t — > — oo) and the initial state of 
the Universe is nearly flat, cold and decoupled: g <§; 1, TZX 2 <C 1. 

The tensorial perturbations in the PBB model satisfy an equation 
similar to Eq. 1)18(1 but with the scale factor a replaced by the pump field 
a = e - ^/ 2 a, that is glj 



Note that a coincides with the scale factor in the Einstein frame. [In the case 
ofD = d + n+ l dimensions with n internal dimensions contracting and 
d external dimensions expanding, isotropically, if tensorial perturbations 
depend only on the external dimensions, then Eq. ((42(1 is still valid but 
the pump field a = e~^l 2 a( d_1 )/ 2 b n l 2 , b being the scale factor of internal 
dimensions .44, •] Let us consider the simplest scenario (henceforth minimal 
model) characterized by an isotropic superinflationary phase with static 
internal dimensions, followed by RD era. In this case, the solution of the 
background field equations is 0Jj a(rj) = (— 7y)( 1_v/ 3)/ 2 and ip(rj) — tpo — 
V3 log(— if), and during RD era, a(r)) = (r) — 2r)*)/(H s r)*) and ip(jj) = ip*, 
where H s ~ M s = go Mpi, with M s the fundamental string mass and go 
the string coupling at present time. Since during the superinflationary era 
the Hubble parameter is not constant (H ^0), in fact it increases in time, 
we expect the GW spectrum for fluctuations re-entering the Hubble radius 
during RD era be non flat. Indeed, using Eq. I(39|) with 71 — 1/2 and 72 = 1, 
it can be easily derived that the spectral slope is (441 1521 15*3*] tit = 3. If 
we assume that CMB photons we observe today carry the (red-shifted) 
energy of the primordial perturbations, i.e. fi 7 (to) = (H*/Hq) 2 (a*/ao) 4 > 
and at the end of the superinflationary phase ~ H s and <?* ~ go, 
then the frequency at the end point of the spectrum is /* = k*/(2irao) — 
a, H s /(2ira ) - 10 10 Hz, and [*fl l5**"l IB***] 



Depending on the present value of the string coupling [""J g = 0.03-0.3, 
the GW spectrum amplitude at the end point frequency is ~ 10 _7 ~10 -5 
[0 7 (*o) ~ 10 -4 ], thus rather close to the BBN bound discussed in Sec- 
tion""""""""""] In the right panel of Fig. [H we show two examples of the stochas- 
tic GW background in the range of frequencies of space- and ground-based 
detectors. In the minimal model the spectrum has the peak in the GHz 
region and being elsewhere rather steep, it cannot be detected by current 



hl( V ) + 2^hl(r ] ) + k 2 h i :(r,)=Q. 



(42) 




(43) 
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Figure 4. In the left panel we sketch the evolution of H and of some physical wave- 
lengths during the superinfiationary PBB phase, RD and MD eras. The shaded box 
refers to the spacetime region around the would-be big-bang singularity for which we do 
not have a complete description, yet. In the right panel we show two examples of the 
stochastic GW spectrum. 

and planned GW detectors. Refs. E31 EOI showed that if the PBB su- 
perinfiationary phase is followed by other cosmological eras where (partial) 
higher-order derivatives and/or quantum loops corrections in the effective 
action <|41|> are included, the GW spectrum can have branches where it 
is less steeper, even flat, and could have amplitudes detectable by space- 
and ground-based detectors. However, robust predictions can be made only 
when a complete description of the transition from pre- to the post-big- 
bang era will be available, including the understanding of the big-bang sin- 
gularity. Indeed, the spectral slope for / > /b in Fig. depends crucially 
on the cosmological backgounds around the would-be big-bang singularity. 

4.2.2. Bouncing- Universe scenarios 

New cosmologies motivated by string theory and influenced also by brane- 
world ideas have been proposed 0H1 EH] ■ As in the PBB model, in those 
scenarios it is assumed that the Universe had a long existence prior to the 
big bang singularity. Here, as an example, we shall focus on the bouncing- 
Universe model discussed in Ref. 0Hj • This model uses two four-dimensional 
boundary branes, one of which is the visible brane where our Universe lives, 
the other is the hidden brane. The two branes are originally separated by 
a finite distance which can be regarded, according to the M-theory conjec- 
ture jnS] , as the dilaton field in the strong-coupling regime. This scenario is 
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built on the following four-dimensional effective action on the visible brane 
(in Einstein frame): 



lD7rGr 



(44) 



The hidden boundary brane (or bulk brane) is assumed to be initially in a 
state with minimal energy, and it is flat, parallel to the visible brane and at 
rest. Non-perturbative effects, for example due to a potential of the kind 
V(4>) = —Vq e~ c< ^ with c > 0, can generate interactions between the visible 
and bulk brane, causing the bulk brane to move toward the visible one and 
collide with it. In this way the kinetic energy of the bulk brane can be 
converted into radiation, and the hot big-bang era starts. 

Differently from the PBB scenario, in the model of Ref. the Uni- 
verse evolves from strong (g 2 = ^> 1) to weak coupling (g 2 = <C 1), 
and both in the string- and Einstein-frame it undergoes an accelerated 
contraction. Moreover, it is assumed that going toward the singularity 
(jj — > 0~) the potential V changes its shape and goes to zero. Thus the cos- 
mological background in the Einstein frame can be roughly described by a 
phase with <2e(?7) ~ (^vY with < e -C 1 (very slow contraction), followed 
by a phase with (almost) zero potential, i.e. a superinflationary phase of 
the PBB kind, with a^(rj) ~ (— if) 1 ! 2 and then a RD era. The pump field 
a entering the equation of tensorial perturbations coincides with the scale 
factor in the Einstein frame (a = tie), so for fluctuations which go out of the 
Hubble radius during the phase dominated by the potential and re-enter the 
Hubble radius during RD era, the discussion around Eqs. gives: 
7i = e > and 72 = 1, thus nx ~ 2 + 2e, whereas for fluctuations that go 
out of the Hubble radius during the phase with (almost) zero potential and 
re-enter during RD era, we recover the PBB result Ut = 3. 



Hence, both the PBB and bouncing-Universe scenarios predict very 
steep GW spectra at very low frequency, so no contribution of tensorial 
fluctuations to the CMBR inhomogeneities. Thus, if a tensorial component 
will be found in CMBR, then the current version of those scenarios should 
be rejected. 

Due to the collapse phase, in both scenarios discussed in this section the 
initial (classical) tensor fluctuations are not de- amplified [HJi as occurs in 
potential-driven inflation (see discussion at the beginning of Section l4.1T [1 . 
This result, though paradoxical, does not imply that the initial value of 
tensor inhomogeneities must be fine tuned to unnaturally small value. In- 
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deed, it can be derived [31] that the energy density of tensor waves is 
indeed de-amplified. However, to solve the homogeneity problem in those 
superstring-inspired models, we have to apply constraints more severe |31] 
than in potential-driven inflation models. 

It is worth to mention that some aspects of both the PBB and bouncing- 
Universe scenarios, have been questioned. Doubts on the naturalness of the 
initial conditions and on the actual solution of the homogeneity and flatness 
conundra were raised |5f>| . The formation of large-scale structures and 
CMBR inhomogeneities cannot be explained in those models as naturally 
as in the standard inflationary models (see Section 13.1(1 using adiabatic 
perturbations. [In the case of the PBB scenario, a way out has been recently 
suggested by converting isocurvature into adiabatic fluctuations [SZ] during 
the post-big-bang phase.] More importantly, a debate [SHj on whether the 
scalar perturbations can be unambiguously calculated in those superstring- 
motivated models is still underway. Those issues can be clarified only when 
a complete and unique description of the transition from pre to post era 
will be available. 

4.3. Non-standard equations of state in post big-bang eras 

In this section we want to give another example of non-flat relic GW spec- 
trum due to the presence, soon after inflation, of a phase which differs from 
the RD era. The possibility of having an all variety of spectra by intro- 
ducing non-standard, i.e. different from RD or MD, equations of state in 
post-big-bang eras was originally pointed out by Grishchuk 27 . 

An interesting example is the case in which the inflationary phase is 
followed by an expanding phase driven by an effective source whose equation 
of state is stiffer than radiation. Peebles and Vilenkin discussed a 
model where the occurrence of a stiff post-inflationary phase is dynamically 
realized through a potential of the kind V(4>) — A (0 4 + M 4 ) for <fi < 
and V((f>) = AM 8 /(</> 4 + M 4 ) for </> > 0, with A - 1CT 14 and M ~ 
10 6 GeV. This model has been denoted quintessential inflation and was 
originally introduced to explain the dark-energy conundrum, that is the 
present accelerated expansion of the Universe. The inflaton field starts its 
evolution at <fi <C — Mpi and rolls toward zero. Inflation ends at <j) ~ —Mp\ 
when significant part of the potential energy has turned into the kinetic 
energy of <f>. The cosmological evolution is then driven by the kinetic energy 
of (f> (the equation of state is p = p), thus the scale factor a(rj) ~ ^Jrj. 
Following the discussion around Eqs. (|3H| - P^|l . we have 72 = 1/2, and 
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Figure 5. We show an example of the stochastic GW spectrum originated in quintessen- 
tial inflation and contrast it with the sensitivity of space- and (correlated) ground-based 
GW detectors. For a detail discussion of the GW spectra as function of the free param- 
eters see Ref. |60| . 



assuming that during quintessential inflation |60| a(rj) ~ — 1/77, we have 
71 = —1. So, the spectral slope for fluctuations that exit the Hubble radius 
during quintessential inflation and re-enter the Hubble radius during the 
stiff era is nr = 1 |591 IK0| . The GW spectrum increases linearly as a 
function of the frequency. An example of such spectrum is shown in Fig. [5] 
where the flat branch refers to gravitons that go out of the Hubble radius 
during inflation and re-enter the Hubble radius during RD era. The peak 
of the spectrum is firmly localized at 100 GHz, with amplitude < 10~ 6 ; the 
frequency fb at which the flat branch starts depend on the free parameters 
111 the model ^ . A very detailed analysis which includes the discussion of 
the free parameters entering the spectrum and the comparison with current 
and planned GW detectors can be found in Ref. 60 . Unfortunately, in 
Peebles and Vilenkin quintessential model |59| the spectrum is few order of 
magnitudes below the sensitivity achievable with advanced ground-based 
detectors. As shown in Ref. [HI], by relaxing the assumption the scalar 
field driving the kinetic-energy phase coincide with the inflaton field, it is 
possible to shift the spectrum ending frequency at lower frequencies and 
increase by few orders of magnitudes the amplitude of the GW spectrum 
in the frequency range of ground-based detectors. 
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5. Gravitational waves from cosmic strings 

Topological defects could have been generated during symmetry-breaking 
phase transitions in the early Universe. Since the beginning of the 80s 
they received significant attention [^2] as possible candidates for seeding 
structure formation. Recently, more accurate observations of CMBR in- 
homogeneities on smaller angular scales and compatibility with density- 
fluctuation spectrum on scales of 100 1 Mpc, restrict the contribution of 
topological defects to at most ~ 10% [S3]. In the following, we shall restrict 
our discussion to cosmic strings formed when a £7(1) local gauge symmetry 
is spontaneously broken [see Ref. 64 for the possibility of producing cosmic 
strings from brane collision at the end of brane inflation scenario 65 ]. 

Cosmic strings are characterized by a single dimensional scale: the mass- 
per-unit-length fi. Since they do not have any ends, they are in the form 
of loops, smaller or larger than the Hubble length. [The length of a string 
is defined as the energy of the loop divided by fj,.] Those loops have very 
large mass-per- unit-length, on the order of 10 22 g/cm if the scale at which 
the symmetry is broken is the GUT scale. Their tension is equal to their 
mass-per-unit-length, so they oscillate relativistically emitting GWs and 
shrinking in size. It is generally assumed, though difficult to prove, that 
a network of cosmic strings did form at some time during the evolution of 
the Universe. In this network the only relevant scale is the Hubble length. 
Small loops (smaller than Hubble radius length) oscillate, emit GWs and 
disappear, but they are all the time replaced by small loops broken off 
very long loops (longer than Hubble radius length). The wavelength of 
the GW is determined by the length of the loop, and since in the network 
there are loops of all sizes, the GW spectrum is (almost) flat in a large 
frequency band, extending from / ~ 10~ 8 Hz to / ~ 10 10 Hz. The GW 
spectrum is nearly Gaussian and has been evaluated in detail in Refs. |66j . 
If r is the characteristic radius of the loop, the quadrupole moment is 
Q ~ \i r 3 and by assuming an oscillation period of r ~ r, the energy 

...2 

emitted is d£ jdt ~ P ~ Gn Q ~ 7 Gn A* 2 1 where 7 is the dimensionless 
radiation efficiency ~ 60. Numerical simulations of cosmic-string network 
would predict ft GW ~ P/p c < 10 _9 -10~ 8 for cosmic strings with GnM < 
10~ 6 . Thus, the Gaussian GW background from cosmic strings could be 
detectable by second and third generation of ground-based interferometers, 
and by LISA as well, being above the galactic stochastic GW background 
of WD binaries. 

More recently, Damour and Vilenkin |67j found that strong beams of 
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high-frequency GWs could be produced at cusps (where the string reaches 
a speed very close to light velocity) and at kinks along the string loop. 
[In time domain they are "spikes" or "bursts".] As a consequence of 
these beams the GW background discussed above becomes strongly non- 
Gaussian. The most interesting feature of these GW bursts is that they 
could be detectable for a large range of values of G^fJ-, larger than the usu- 
ally considered search for the Gaussian spectrum. For example, if the aver- 
age number of cusps in a string oscillation is only 10%, the GW bursts could 
be detectable by ground- and space-based detectors up to GnA* ~ 10 -13 . 
If the number of cusps in a string oscillation is very small, kinks, which 
have smaller bursts but are ubiquitous, could generate a signal detectable 
by LISA for a large range of values of G^fi [see Figs. 1, 2 in Ref. |H7)]. 

The most stringent bound on the GW backround from cosmic strings 
comes from pulsar timing [see Section I3.4.3j . The analysis of Ref. |o"7] 
suggests that the GUT value GnA«gut ~ is still compatible with 

existing pulsar data, and if accuracy will be improved, pulsar timing could 
allow to detect the GW bursts up to GnA* ~ 10 -11 . 

6. Gravitational waves from other mechanisms occurring in 
early Universe 

During its history the Universe could have undergone phase transitions, 
corresponding to symmetry breaking of particle-physics fundamental inter- 
actions. 

In a first-order phase transition the Universe is initially trapped in a 
high-temperature metastable phase (unbroken-symmetry phase) , where the 
false vacuum is separated from the true vacuum by a barrier in the potential 
V(4>, T) of some scalar field (f> driving the transition {(f) is also called the 
order parameter) [see Fig HJ . The transition from the metastable to the 
ground state takes place via quantum tunneling across the barrier with 
random nucleation of bubbles: true vacuum is inside the bubbles, false 
vacuum is outside them. If the temperature when bubbles form is too high, 
they are born small, the volume energy cannot overcome the shrinking 
effect of the surface tension and they disappear quickly; but when the 
temperature drops below a critical temperature, bubbles are born larger, 
and since the latent heat released during the transition is converted to 
bubble-wall kinetic energy, they expand, approach velocity close to the 
speed of light, collide and leave the Universe in the broken-symmetry phase. 
When the collision occurs, the spherical symmetry is broken and GWs, as 
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Figure 6. We show the typical temperature-dependence of a potential V(tf>,T) in a 
first-order phase transition <j> being the order parameter. At the temperature T the true 
(< <j> >^ 0) and false (conventionally < <fc >= 0) vacua are degenerate, and at T the 
false vacuum becomes unstable. The transition occurs at some temperature T < T + < T 



well as other particles are radiated away Since the Universe is 
expanding, the temperature T* at which the transition takes place can be 
obtained comparing the probability of bubble nucleation per unit time and 
volume with the expansion rate of the Universe at that temperature. The 
transition occurs when the probability for a single bubble to be nucleated 
within one horizon volume is on the order of one. 

Two parameters determine the stochastic GW spectrum [JO] : the bubble 
nucleation rate per unit volume (3 (Y — To e' 3 ') and the ratio a between the 
false vacuum energy density and the energy density of the radiation at the 
transition temperature T*. The bubble collision produces a GW spectrum 
which is strongly peaked at the frequency characteristic of the nucleation 
rate, i.e. 27r/ pca k ~ (3. A detailed analysis gives |7H] : 

/P .^5.2X10-»(|-)( I ^-)(^) 1 %,, (45) 

where and <?* are the Hubble parameter and the number of degrees of 
freedom at the time of transition. Typical values for an electro- weak phase 
transition (EWPT) are 0/H* ~ 10 2 -10 3 , = 0(100) GeV, for which 
/peak — 10~ 4 -5 x 10~ 3 Hz, and lies in the LISA frequency band. The GW 
spectrum was calculated in Ref. [70]; it reads 

where k quantifies the fraction of latent heat that is transformed into 
bubble-wall kinetic energy and Vb is the bubble expansion velocity. At 
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frequencies lower than / pea k the GW spectrum cx f 2 S [70], while at higher 
frequencies it drops off as cx J -1 ' 8 |70] . 

Non-perturbative calculations done using lattice field theory have shown 
that there is no first-order EWPT in the SM of particle physics for Higgs 
mass larger than W masses (current results predict an Higgs mass larger 
than W masses). In Minimal Supersymmetric Standard Models (MSSM), if 
the Higgs mass is in the range 110-115 GeV, there is the possibility of having 
a first-order phase transition, if the right-handed stop mass has a mass 
in the range 105-165 GeV [the stop particle is the scalar supersymmetric 
partner of the top quark] . However, in this case the GWs produced are too 
weak. For example, for a Higgs mass of 110 GeV and right-handed stop 
mass of 140 GeV, there exist regions in the parameter space where |7J 
a ~ 10~ 2 , k ~ 0.05 and H*/fi ~ 10~ 4 and thus h% il GW ~ 10~ 19 around 10 
mHz. The strength of the transition can be enhanced in next-to-minimal 
supersymmetric standard models (NMSSM). Apreda et al. |7J investigated 
the NMSSM obtained adding a gauge singlet in the Higgs sector ,72 . This 
model is rather attractive also because it can explain the observed baryon 
asymmetry. Apreda et al. found that there exist regions in the parameter 
space for which the GW spectrum is h$ Hqw ~ 10 _15 -10 -10 at / pC ak — 
10 mHz. Note that for frequencies 10 _4 -3 x 10 -3 Hz the stochastic GW 
background from WD binaries could "cover" the GW spectrum from first- 
order phase transitions. If this GW background will be ever detected it will 
be a signature of supersymmetry and, combined with experimental bounds 
from future particle colliders, it could allow to discriminate between various 
supersymmetric models. 

A stochastic GW background could be also produced during a first- 
order phase transition from turbulent (anisotropic) eddies generated in the 
background fluid during the fast expansion and collision of the true- vacuum 
bubbles [231711 EH In the NMSSM [7J] there exist regions of the param- 
eter space where [7J /igSlcw ~ 10~ 10 with peak frequency in the mHz. 
Recently, the authors of Ref. [7J] evaluated the stochastic GW background 
generated by cosmic turbulence before neutrino decoupling, i.e. much later 
than EWPT, and at the end of a first-order phase transition if magnetic 
fields also affect the turbulent energy spectrum. The observational perspec- 
tives of those scenarios are promising for LISA. 

Turner and Wilczek |75| pointed out that if inflation ends with bub- 
ble collisions, as in extended inflation, the GW spectrum produced has a 
peak in the frequency range of ground-based detectors. Subsequent analy- 
ses have shown that in two-field inflationary models where a field performs 
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the first-order phase transition and a second field provides the inflationary 
slow rolling (so-called first-order or false vacuum inflation |7fi|). if the phase 
transition occurs well before the end of inflation, a GW spectrum peaked 
around 10-10 3 Hz, can be produced |77j . with an amplitude large enough, 
depending on the number of e-foldings left after the phase transition, to be 
detectable by ground-based interferometers. A successful detection of such 
a spectrum will allow to distinguish between inflation and other cosmologi- 
cal phase transitions, like QCD or electroweak, which have a different peak 
frequency. 

Another mechanism that could have produced GWs in the early uni- 
verse is parametric amplification after preheating T8: . During this phase 
classical fluctuations produced by the oscillations of the inflaton field </> can 
interact back, via parametric resonance, on the oscillating background pro- 
ducing GWs. In the model where the inflaton potential contains also the 
interaction term ~ 2 X 2 j X being a scalar field, the authors of Ref. |75) 
estimated f^cw ~ 1CP 12 at / m i„ ~ 10 5 Hz, while in pure chaotic inflation 
S1 GW < 1CT 11 at / min - 10 4 Hz. [See Fig. 3 in Ref. [75] for the GW 
spectrum in the range 10 6 -10 8 Hz.] Unfortunately, the predictions lie in 
frequency range where no GW detectors have been planned so far, although 
EM cavities to detect GWs were proposed |T5] . 

7. Gravitational waves in brane world scenarios 

During the last years there have been feverish activities around brane- world 
models [SU||^. Those scenarios are based on the idea that large (> l/Mpj) 
spatial, extra dimensions might exist. They assume that ordinary matter 
is confined onto a three-dimensional subspace, called the (visible) brane, 
which is embedded in a larger space, called the bulk. Only gravitational 
interactions are allowed to propagate in the bulk. If n is the number of 
extra dimensions, then the four-dimensional Planck mass Mpj is a derived 
quantity, while the fundamental scale is determined by the gravitational 
mass Mfund in n + 4 dimensions. In the case of a flat bulk, / being the 
typical length of the extra dimensions, M™ = M f ™+ 2 l n with Mf un d ranging 
between TeV to Mp\. In the following, we shall mention results obtained 
in five-dimensional brane- world models [n = 1). 

By normalizing properly Eq. 132(1 and taking the limit kij <C 1, it is 
easily derived that long-wavelength GWs generated in de Sitter inflation 
have a (almost) constant amplitude fc 3 / 2 hk = H^/Mvi, which depends 
on the two scales in the problem: the Hubble parameter during inflation 
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and the low-energy Planck mass. In five-dimensional brane-world models 
there is an additional scale in the problem, the scale I associated to the 
fifth dimension. Quite generically, we could expect [82) that if inflation on 
the (visible) brane occurs at scales smaller than the curvature of the extra 
dimension, i.e. H I 3> 1, then the amplitude of GWs f can be modified with 
respect to the value predicted in standard four-dimensional theory - for 
example [22 M E3 k 3 / 2 h k = f(Hl)H/M Pl with some function /. « 
In the brane-world models analysed in Refs. [EUl 1841 [S2j , it was found that 
for ff I > 1 the tensor amplitude is enhanced with respect to the standard 
result, i.e. f(lH) > 1. Moreover, the ratio between scalar and tensorial 
perturbations T/S may also differ from the value predicted by standard 
four-dimensional theory E2] and it could depend on specific features of 
the brane-inflation model. However, as seen in Section [3 the overall shape 
of the relic GW spectrum depends not only on the evolution of the tensor- 
mode amplitude during inflation but also on the post-big-bang phases, i.e. 
on the era at which the mode re-enters the Hubble radius. Only when 
those cosmological phases will be consistently described in the brane-world 
scenarios, we will have robust predictions for the GW spectrum. 

By implementing quintessential inflation on the visible brane, and evalu- 
ating the Bogoliubov coefflents all along the Universe evolution, the authors 
of Ref. [SS] found results similar to the ones discussed in Section 14.31 the 
GW spectrum has a branch where it increases linearly as function of the 
frequency. 

GWs can be also produced by excitations of the so-called radion field, 
which controls the size of the extra dimension, and by inhomogeneities in 
the displacement of the brane [HZj • Those waves should peak at a frequency 
fixed by the scale of the extra dimension. LISA could observe excitations 
emitted at scales from millimeters to microns, while ground-based interfer- 
ometers could observe rather small extra dimensions, up to ~ 1CP 12 mm. 

In some brane world scenarios the causal propagation of gravitational 
and luminous signals can be different [HH1 EH] ■ In the case of an asymmetric 
warped spacetime of the kind (EH], 

ds 2 = -n 2 (l) dt 2 + a 2 (I) ( 2 + r 2 dti 2 ) + b 2 {l) dl 2 , (47) 



Here the GW refers to the homogeneous Kaluza-Klein massless mode. Massive Kaluza- 
Klein gravitons decay very fastly 1831 1841 1821 . so their production is very suppressed 
during brane inflation. 

g In the models analysed by Frolov and Kofman 1821 . f(Hl)/Mp\ = 1/Mpi; n fl where 
Mpi i n fl is the Planck mass during inflation. 
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where I is the coordinate along the extra dimension, and k — ±1,0 is the 
spatial curvature of the three-dimensional sections parallel to the brane, 
the local speed of light c(l) = n(l)/a(l) depends on the position I along 
the extra dimension. Since GWs can propagate in the bulk, a GW signal 
emitted at point A on the brane can take [if c(l) is increasing away from 
the brane] a short-cut in the bulk, and to an observer located at point B 
on the brane it appears quicker than a photon traveling along the brane 
from A to B. If true, this effect predicts a difference between gravitational 
and luminous speed. By detecting GWs and EM waves from 7-ray bursts 
and supernovae, ground-based detectors of second generation can put upper 
limits on 8c/ c at most on the order of 10 -17 ^2] for 7-ray bursts, and 10 1 
for a supernova in the Virgo cluster of galaxies. 

8. Extraction of cosmological parameters from detection of 
GWs 

In this section we want to discuss very briefly another possible way GW 
detection can be used to probe cosmology, although not directly the physical 
mechanisms occurring in the very early Universe. 

It was realized long ago that binaries made of compact bodies, like BHs 
or NSs, which spiral in toward each other loosing energy because of the 
emission of GWs, are "standard candles" [5U]. Indeed, once the binary 
masses and spins, position and orientation angles are specified, the GW 
signal passing-by the detector depends only on the luminosity distance 
dij(z, Jljif, Oa, • • •) between binary and detector. By using three ground- 
based interferometers or LISA (and suitably high signal-to-noise ratio) it 
is possible to determine the location and orientation of the binary in the 
sky, extract the masses, the spins and the cosmological distance but not 
the source's cosmological red-shift. Indeed, the inspiral GW signal emitted 
from a binary with masses (mi, m^) at red-shift z cannot be distinguished, 
except for an overall factor in the amplitude, from the one emitted from a 
local binary with masses [(1 + z) mi, (1 + z) m2\- To break the degeneracy 
and obtain the distance-red-shift curve and the cosmological parameters, 
one could associate the binary-coalescence event to an EM event which 
had very clear emission or absorption lines, from which one could read z. 
[Because the errors in position determinations can be rather large, the ex- 
istence of EM counterpart can also improve the accuracy in measuring the 
luminosity distance - for example by mapping massive BH binaries with 
LISA (Wl/^l decreases from 1-10% to 0.15-1% |5T].] However, Holz and 



Hughes [2] recently pointed out that, practically, because of gravitational- 
lensing effects (GWs are lensed as EM radiation are lensed) the precision 
with which the luminosity distance can be determined is degraded, though 
still slightly better than what Type-la supernovae probes can reach, but 
unless the event rates is high enough the extraction of cosmological param- 
eters cannot be done with good accuracy. 

9. Summary 

As seen, the search for GWs from the very early Universe is rather chal- 
lenging but the outcome is certainly worth the effort. 

Theoretically, so little we know about the very early evolution of our 
Universe, that the predictions for the stochastic GW background from stan- 
dard inflation [37J EH] and/or superstring-motivated models H7JHB]) includ- 
ing also brane- world scenarios [HS1 EU IHH EHl EZj > should be considered 
just as indications. What we learned from those models, is that it could 
well be the GW spectrum is not (almost) flat all the way from / ~ 1CP 16 
Hz to / ~ 10 10 Hz. There could be frequency regions in which it increases 
or decreases [27] due to cosmological phases existed before the would-be 
big-bang singularity, and/or post-big-bang phases with equations of state 
different than radiation or matter ones, whose presence we cannot currently 
exclude. [Even an initial state different from vacuum could affect the relic 
GW spectrum |53|. as seen from Eq. (|T7Jl .] Future CMB polarization experi- 
ments | 43| could detect the tensor contribution at very large wavelength and 
discriminate within plethora of inflationary models. The detection of GWs 
from bubble collision EH3 EH E3 C7J and/or turbulent motion [Ml fTTl 1751 
174) in the primordial plasma, and especially the frequency at which the GW 
spectrum is peaked, will reveal if first-order phase transitions occured, and 
at which particle physics mechanism they are associated to: QCD, EW, last 
stages of inflation, etc. If a network of cosmic strings ever formed, it should 
have produced a Gaussian [SB] or even strongly non-Gaussian |B7j GW spec- 
trum, which could be detectable by the second generation of ground-based 
interferometers and LISA. If not observed, we could put upper limits on 
few free parameters. An independent estimation of the cosmological pa- 
rameters could be obtained by detecting GWs from coalescing binaries at 
fairly high red-shift |SDl EU ■ 

Experimentally, since the GW signal from astrophysical sources is much 
more promising than the one from the early Universe, the ground- and 
space-based detectors were conceived and designed to detect the formers. 
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Various scientists in the GW community have already started thinking at 
future detectors, like the follow-on LISA missions which could focus more 
on early-time cosmological signals [23 EH1 EHD 013 EI] • The technological 
challenges for these missions are considerable and deserve very careful in- 
vestigations. Finally, ground-based GW detectors in the high-frequency 
region of MHz 79 and GHz [22 have gained more attention j4Tl ITS] . 

So, even if the detection of GWs from primordial Universe is still some- 
what far ahead of us, maybe, if nature is kind with us, we will not wait for 
a long time. There could be surprises. 
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